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Abstract 

The impact of the leading quantum gravity effects on the dynamics of the Hawking evaporation 
process of a black hole is investigated. Its spacetime structure is described by a renormalization 
group improved Vaidya metric. Its event horizon, apparent horizon, and timelike limit surface 
are obtained taking the scale dependence of Newton's constant into account. The emergence of a 
quantum ergosphere is discussed. The final state of the evaporation process is a cold, Planck size 
remnant. 
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I. INTRODUCTION 



One of the very remarkable features of black hole radiance [l| is the observation that the 



0is 



)y the evaporation process is far more 
. Even a Schwarzschild black hole, 



global spacetime structure of a black hole losing mass 
complicated than that of its static counterpart P, 0, 
when it radiates, does not have a single horizon that fully characterizes its structure, and 
one must distinguish at least three important horizon-like loci. The future event horizon 
(EH) is the boundary of the causal past of future null infinity, and it represents the locus 
of outgoing future-directed null geodesic rays that never manage to reach arbitrarily large 
distances from the hole. The apparent horizon (AH) is defined as the outermost marginally 
trapped surface for the outgoing photons. Classically it can be null or spacelike, in presence 
of quantum radiance it can be timelike also, when regarded as a 3-dimensional surface. The 
third important locus is the timelike limit surface (TLS) or "quasi-static limit" which is 
defined as the locus where static observers become lightlike. The TLS can be null, timelike, 
or spacelike [2]. For a classical Schwarzschild black hole (which does not radiate), the three 
surfaces EH, AH, and TLS are all identical. Upon "switching on" the Hawking evaporation 
this degeneracy is partially lifted. According to the analysis by York P, ^ the AH continues 
to coincide with the TLS for a spherically symmetric emission, but the EH is different from 
AH=TLS. 

In particular, if we approximate the stress-energy tensor near the horizon as a radial 
influx of negative energy which balances the outward Hawking flux at infinity, the event 
horizon is located inside the AH |5j, the portion of spacetime between the two surfaces 
forming the so-called "quantum ergosphere". This name stems from the analogy with the 
classical (stationary) Kerr black hole for which EH=AH7^TLS. Here the ergosphere is the 
space between "the" horizon EH=AH and the TLS, usually called the "static limit". In both 
cases particles and light signals can escape from within the ergosphere and reach infinity. 

The definition of the EH via the locus of outgoing photons that can never reach large 
distances from the hole has the unfortunate "teleological" property of requiring knowledge 
of the entire future history of the hole P, 0| . In particular when the black hole radiance is de- 
scribed semiclassically (quantized matter in a classical geometry), the Bekenstein-Hawking 
temperature and the luminosity diverge for M — > 0, as Tbh oc 1/M and L oc re- 
spectively. As a result, this approximation breaks down for very light holes, and in order to 
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determine the final state of tlie evaporation process a mucli more precise treatment, including 
backreaction and quantum gravitational effects, is required. 

In York's work |^ |^ , which is strictly within the semiclassical approximation, the "tele- 
ological" problem is circumvented by relaxing the definition of the EH in the following way. 
Rather than demanding that the photons "never" reach infinity he demands only that they 
are imprisoned by the event horizon for times which are very long compared to the dynamical 
time scale of the hole. Using this working definition of the EH he is then able to determine 
its location to first order in the luminosity L. In this manner the difficult question about 
the real final state of the evaporation is not touched upon. 

It is the purpose of the present paper to analyze the dynamical evaporation process and 
the corresponding spacetime structure of a radiating Schwarzschild black hole. We include 
the leading quantum gravitational corrections of the geometry which, as we shall discuss, 
seem to lead to a termination of the evaporation process and the formation of a cold, Planck 
size remnant. Our main tool will be the "renormalization group improvement" of classical 
solutions, a technique which is very popular in conventional field theory. 

In fact, recently a lot of work went into the investigation of the nonperturbative 



renormalization group (RG) behavior of Quantum 
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26| and its possible manifestations 
28[, a "RG-improvement" 



37l |. In particular, in 

of the Schwarzschild metric has been performed and the properties of the corresponding 
"quantum black hole" have been explored. The improvement was based upon the scale de- 
pendent ( "running" ) Newton constant G{k) obtained from the exact RG equation for gravity 
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421. Here k 



[7| describing the scale dependence of the effective average action 
denotes the mass scale of the infrared cutoff which is built into the effective average action 
rfc[5'^i/] in such a way that it generates the field equations for a metric which has been aver- 
aged over a spacetime volume of linear dimension k~^. The running of G is approximately 
given by 

G{k) = , , ^,1 (1.1) 



1 + Go A;2 

where Gq denotes the laboratory value of Newton's constant, and u; is a constant. At 
large distances {k 0), G{k) approaches Gq, and in the ultraviolet limit {k oo), it 
decreases as G{k) oc l/k"^. This is the fixed point behavior responsible for the conjectured 
nonperturbative renormalizability of Quantum Einstein Gravity 3, 



ii|,Q. 



In the RG improvement scheme of '28l] the information about the fc-dependence of G is 
exploited in the following way. The starting point is the classical Schwarzschild metric (in 
Schwarzschild coordinates) 

ds^ = -f{r)dt^ + f{r)-^dr^ + r^dil^ (1.2) 

with dfl"^ = dO"^ + sin^ 9d(f>'^ and the classical lapse function /(r) = 1 — 2GoM/r = fdassi^). 
The RG improvement is effected by substituting, in /ciass('"), Gq by the r-dependent Newton 
constant G{r) = G{k = k{r)) which obtains from G{k) via an appropriate "cutoff identifica- 
tion" k = k{r). In flat space the natural choice would be oc 1/r. In we argued that in 
the Schwarzschild background the correct choice, in leading order at least, is k{r) = ^/d{r) 
where ^ is a constant of the order of unity, and d(r) = /q'^ (ir'|/ciass('"')l^^^^ is the pro per 
distance from a point with coordinate r to the center of the black hole. (We refer to [28| 
for a detailed physical justification of this choice.) While the integral defining d{r) can be 
evaluated exactly, it is sufficient to use the following approximation which becomes exact 
for both r ^ oo and r — > 0: 

The resulting G{r) = G{k = ^/d{r)) reads 

G(r) = ^ (1 4) 

r3 + uGo[r + ^GoM] ^ ^ 

where u = ooC,'^. In these equations the parameter 7 has the value 7 = 9/2 if one sets 
k = ^/d{r) as above. It turns out, however, that most of the qualitative properties of the 
improved metric, in particular all those related to the structure of its horizons, are fairly 
insensitive to the precise value of 7. In particular, 7 = (corresponding to k = C,/r) and 
7 = 9/2 where found [28] to lead to rather similar results throughout. For this reason 
we shall adopt the choice 7 = in the present paper. It has the advantage that with this 
choice many calculations can be performed analytically which require a numerical treatment 
otherwise. 

The metric of the RG improved Schwarzschild black hole is given by the line element 
with 

(1.5) 
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Let us briefly list its essential features^, 
a) There exists a critical mass value^ 

McT = ^/uj/Go = mpi (1.6) 

such that /(r) has two simple zeros at r_ and r+ > r_ if M > M^r, one double zero at 



r_ = y/ooGo if M = M^r, and no zero at all if M < M^. For M > the zeros are at 



r± = GoM [1 ± v^r^] (1.7) 

with the convenient abbreviation 

n^^ = u(^Y (1.8) 

The spacetime has an outer horizon at r+ and in inner (Cauchy) horizon at r_ . At Mcr, 
the black hole is extremal, the two horizons coincide, and the spacetime is free from any 
horizon if the mass is sufficiently small, M < M^- 

b) The Bekenstein-Hawking temperature Tbh = nf^ir is given by the surface gravity at 
the outer horizon, k = |/'(r+). Explicitly, 

T (M) ^ ^^^^ ^ ^/n{l-n) M,, ^nii-n) 

''''^ ' AtxGqM 1 + v^r^ 47rGoMer 1 + ^TT^ Atiu 1 + ^/l^ ^ ' ^ 
This temperature vanishes for M \ M^, i.e. Q y 1, thus motivating the interpretation 
of the improved Schwarzschild metric with M = Mcr as describing a "cold" remnant of the 
evaporation process. 

c) The energy flux from the black hole, its luminosity L, can be estimated using Stefan's 
law. It is given by L = aA{M)T^i{{M)'^ where a is a constant and A = Anr"^ denotes the 
area of the outer horizon. With ()1.7|) and ()1.9|) we obtain 

^ ^ (47r)3cu2 [i + ^/r^]2 ^ ^ 

For a single massless fleld with two degrees of freedom one has a = n'^/QO. (We use units 
such that h = c = kB = I.) 



^ All formulas quoted refer to 7 = 0, but the qualitative features are the same for 7 = 9/2; see |23| for 
details. 

^ We define the (standard) Planck mass and length in terms of the laboratory value Gq: mpi = ^pj^ = 



II. THE QUANTUM-CORRECTED VAIDYA METRIC 



Our aim is to find a metric which describes the history of an evaporating Schwarzschild 
black hole and its gravitational field. In the small luminosity limit {L 0) this metric is 
supposed to reduce to the static metric of the RG improved Schwarzschild spacetime. 

We begin by reexpressing the metric ()1.2j) with the improved lapse function ()1.5|) in terms 
of ingoing Eddington-Finkelstein coordinates {v,r,6,(j)). We trade the Schwarzschild time t 
for the advanced time coordinate 

v = t + r*, = I dr'/f{r') (2.1) 



Here r* is a generalization of the familiar "tortoise" radial coordinate to which it reduces if 
G{r) = const. For G{r) ^ const the function r* = r*(r) is more complicated, but its explicit 
form will not be needed here. Eq. (|2.1|) implies dv = dt + dr/ f{r), turning (|1.2p with (|1.5|) 
into 

ds"^ = -[1 - 2G{r)M/r] dv'^ + 2dvdr + r^dfi^ (2.2) 

Eq. ()2.2p is exactly the Schwarzschild metric in Eddington-Finkelstein coordinates, with Gq 
replaced by G{r). It is thus reassuring to see that the two operations, the RG improvement 
Go G{r) and the change of the coordinate system, can be performed in either order, they 
"commute" . 

28[ and summarized in the previous section 



The thermodynamical properties derived in 
refer to the metric ()2.2|) . In the exterior of the hole the spacetime is static, and while we 
can deduce a temperature and a corresponding luminosity from its periodicity in imaginary 
time (or by computing the surface gravity at r+ directly) the backreaction of the mass-loss 
due to the evaporation is not described by (j2.2|) . From the static metric we obtained the 
mass dependence of the luminosity, L = L{M). Using this information we can compute the 
mass of the hole as seen by a distant observer at time v, M{v), by solving the differential 
equation 

—^M(v) = LiMiv)) (2.3) 
dv 

In our case L{M) is given by Eq. (jl.l(J|) . To first order in the luminosity, the metric which 
incorporates the effect of the decreasing mass is obtained by replacing the constant M in 
(123) with the M{v) obtained from Ea.dO): 

ds^ = -[1 - 2G{r)M{v)/r] dv^ + 2dvdr + r^dQ^ (2.4) 



For G{r) = const, Eq. ()2.4|) is the Vaidya metric which frequently had been used to explore 



43|, 



4^ |45l. It is a solution of 



the influence of the Hawking radiation on the geometry 
Einstein's equation G^,^ = SttGqT^i, where T^^, describes an inward moving null fluid. In this 
picture the decrease of M is due to the inflow of negative e nerg y, as it is appropriate if the 
field whose quanta are radiated off is in the Unruh vacuum 46l |. 

The metric ()2.4|) can be regarded as a RG improved Vaidya metric. It encapsulates two 
different mechanisms whose combined effect can be studied here: the black hole radiance, 
and the modifications of the spacetime structure due to the quantum gravity effects, the 
running of G in particular. 

It is instructive to ask which energy-momentum tensor T^,^ would give rise to the improved 
Vaidya metric (j2.4p according to the classical equation G^'^ = StxGqT^" . Computing the 
Einstein tensor of ()2.4j) one finds that its only non-zero components are 

r«..rv^-«l (2.5a) 

STrGor^ 



Here the prime (dot) denotes a derivative with respect to r{v). The non-zero components 
()2.5|) contain either r- or f-derivatives but no mixed terms. The terms with r-derivatives of 
G, also present for M{v) = const, describe the vacuum energy density and pressure of the 
improved Schwarzschild spacetime in absence of radiation effects. They have been discussed 
in 3] already. Allowing for M{v) ^ const, the new feature is a nonzero component T**^ ^ 
which, for M < 0, describes the inflow of negative energy into the black hole. 

Taking advantage of the luminosity function L{M), Eq. p.l()|l . we can solve the differ- 
ential equation ()2.3|) numerically and obtain the mass function M = M{v). (We have set 
cr/(4'7r)^u5 = 1 in the numerical calculations in order to reach the almost complete evapora- 
tion for V ~ 200 in units of Vcr-) The result is shown in Fig.(^ for various initial masses, in 
the domain f > 0. In fact, for definiteness we assume that the black hole is formed at f = 
by the implosion of a spherical null shell j^]. Hence M{v) is given by Fig.([Q) together with 
M = for f < 0. We observe that, for any initial mass, M{v) approaches the critical mass 
Mcr for V oo. This behavior is the most important manifestation of the quantum gravity 
effects: according to Eq. ()1.9|) . the temperature Tbh(^) goes to zero when M approaches 
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FIG. 1: The ratio M/Mcr as a function of f/rcr for various initial masses. 

Mcr from above. Hence the luminosity vanishes, too, the evaporation process stops, and 
M{v) ~ Mcr remains approximately constant at very late times, v ^ M~^. 

In Fig.Q and Fig.© we plot the advanced time dependence of the temperature T-^niy) = 
TBH(M(f )) and the luminosity L{v) = L{M{v)), respectively. They are obtained by inserting 
the numerical solution of Eq. ()2.H|l into and 

Both the very early and the very late stages of the evaporation process can be described 
analytically. A black hole with M{v = 0) ^ starts in what we call the "Hawking 
regime" . It is defined by the approximation i7 ^ which is realized if the hole is very heavy 
(M ^ Mcr) or if a; = in the semiclassical limit where the quantum gravity corrections are 
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FIG. 2: Time dependence of the Bekenstein-Hawking temperature during the evaporation process 
for the same initial masses as in Fig.lfH). 



"switched off" . In the Hawking regime, ()1.9p and ()1.10p reproduce the familiar results 

1 



Tbh(M) 
L{M) = 



SttGo M 
B 



B 



a 



(2.6a) 
(2.6b) 



4(47r)3 

It is easy to solve the differential equation —M = L{M) for the luminosity ()2.6b|) . With the 
initial condition M{v = 0) = Mq the solution reads 

1/3 



M{v) = - 3(5/^^) V 



(2.7) 



This is the mass function during the early stages of the evaporation process, valid as long as 
M{v) is well above the critical mass. If one naively extrapolates ()2.7|) to small masses one 
finds M(fo) = 0, implying a final "explosion" with T oo and L oo, after a finite time 
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FIG. 3: The black hole's luminosity as a function of w/rcr, for the same initial masses as in Fig.([T|l 
and Fig.©. 

fo = GoM?/(3-B)- As a consequence of the quantum gravity effects, this is not what really 
happens, however. 

The final part of the evaporation process, where the cold remnant forms, is in the "critical 
regime". It is described by those terms in the above expressions which are dominant for 
M \ Mcr, oi Vt /" 1. From ()1.9|) and p.lO|) we obtain in this approximation: 



Tbh(M) = ^^M^-Ml 



L{M) 



a Go 



(47ru;)3 

Solving -M = L{M) with (|I!8b|l one finds 

M{v) = Mcr + 



- Ml 



(2.8a) 
(2.8b) 



Ml - Mc: 



l + a(Mi - M^,){v-vi) 



(2.9) 
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Here a = cr/(167r^cD^), and vi is a time, already in the critical regime, where M{vi) = Mi 
is imposed. For v oo, the difference M{v) — M^r vanishes proportional to as a 
consequence of which T-Q^iiy) — > and L{v) 0. 

We mentioned already that the RG improved Vaidya metric (j2.4j) can be a correct de- 
scription only to first order in L. In fact, deriving the surface gravity and luminosity from 
()2.4|) the results differ from those for the improved Schwarzschild metric by terms due to the 
v-dependence of M. In our approximation those terms are neglected as they would contain 
additional factors of M = —L. 



III. APPARENT HORIZON AND TIMELIKE LIMIT SURFACE 

Next we turn to the various horizon-like loci of the improved Vaidya metric. Regarded 
as 3-surfaces, all of them are histories of spherical 2-surfaces. 

The apparent horizon is a marginally trapped surface. We determine it from the condition 
that one of the congruences of radial null geodesies, in affine parametrization, has vanishing 
expansion scalar there, = 0. The improved Vaidya metric ()2.4j) has the structure 

ds^ = [-/(r, v)dv + 2dr]dv + r^dn'^ (3.1) 

Along outgoing radial null geodesies we have fdv = 2dr. Hence, parametrizing them as 
r = r(t>), they obey the differential equation f(t>) = /(r(t>), t>)/2 where the dot denotes a 
derivative with respect to v. We can rewrite this equation in the autonomous form 

^x^{X) = u^{x{X)) (3.2) 

with the null vector field 

^ (M^ M^ u^) = (1, 0, 0) (3.3) 

A short calculation reveals that the geodesic equation holds in the form u'^D^u^ = /Cm^ with 
a nonzero function 

/C = ^drf (3.4) 

lence the parameter A in (13.21) is not an affine one. In order for the standard discussion 

in 

6|, |47( to be applicable we must reexpress the null geodesies in terms of an affine parameter 
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A*. Given a solution x(A) of ()3.2p we compute the function A*(A) by integrating 

^X,{X)=expj rfA'/C(x(A')) (3.5) 
and determine its inverse A = A(A,f). Then we define a;^(A^,) = a;^(A(A^,)) which satisfies 

'^-x^(A,) =n^(x*(A,)) (3.6) 



Here 

n''{x) = e-^^'^^u'^ix) (3.7) 
is a new vector field, with T{x) satisfying 

u^'d^T = K. (3.8) 

Using ()3.8p one easily verifies that n^DyU^^ = 0, implying that A* is an affine parameter j^. 

The expansion scalar which determines the location of the AH is the divergence of n'^. 
Using ()3.8|) we find 

e = D^n^" = [Df.u" - /C] . (3.9) 
For f{r,v) = 1 — 2G{r)M{v)/r we have 

DX=-/ + ^^r/ (3.10) 

r 2 

Together with ()3.4p this yields the expansion scalar 

e = K-^f (3.11) 

This is the same result as for the classical Vaidya metric; the r-dependence of G did not 
lead to extra terms. 

Eg. ()3.1H1 tells us that O vanishes if, and only if, / vanishes. According to point a) of 
the Introduction this is the case at r = r_(. and r = r_ with r± defined by ()1.7j) with ()1.8j) . 
(We assume that the r-dependence of G is given by (jl.4j) with 7 = 0). Since r± depends on 
M, it has become a function of the advanced time v now: r±{v) = r±{M{v)). Defined as 
the outermost trapped surface, the AH is characterized by the implicit equation r = r+(f) 
where, explicitly. 



r+{v) 



M{v) , lfM{v)\'i ^ 
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(3.12) 




FIG. 4: Light rays of the outgoing null congruence for M/Mcr = 2. The thick solid line is the 
EH determined numerically, the dashed lines are the outer and inner TLS=AH, and the thin solid 
lines are generic null rays, inside and outside the EH. 

The second horizon surface, the TLS, is defined as the locus where the 4-velocity of 
static observers -u" oc 5" becomes lightlike, with d/dr a vector orthogonal to the r = const 
hypersurfaces. Since we consider a spherically symmetric spacetime, this vector is precisely 
d/dv and AH and TLS coincide in this case, being u^u^gap — Qw — ~f- 

Strictly speaking, apart from the outer AH r = r+{v) there exists also an inner TLS=AH 
r = r-{v), at which the vector field d/dv switches back from spacelike to timehke. 
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IV. THE EVENT HORIZON 



The radial light rays r = r{v) of the outgoing null congruence are to be found by solving 
the differential equation r{v) = f{r{v),v)/2, or explicitly, 



Depending on their initial points, light rays can, or can not, escape to infinity for u — > oo. By 
definition, the "separatrix" separating those two classes of solutions is the event horizon, the 
outermost locus traced by outgoing photons that can never reach arbitrarily large distances. 

It is easy to understand why for a radiating (as opposed to an accreting) black hole the 
EH is inside the TLS. Inside (outside) the TLS, / is negative (positive), implying that the 
light ray's r{v) decreases (increases) if r(t>) < rTLs('y) {r{v) > TTLsiv)). Under certain 
conditions it can happen that due to the hole's mass loss the radius of the TLS decreases 
faster than r{v). As a result, the light ray intersects the TLS, with r = there, and then 
escapes from the hole with r > 0. 

The situation is illustrated in Fig.Q. The region between teh and ttls forms the hole's 
"quantum ergosphere" which owes its existence entirely to the evaporation process. 

Clearly a determination of the EH's VEniv) requires knowledge of M{v) for arbitrarily 
late times v even. Therefore the semiclassical approximation is not sufficient to find the EH 
since it breaks down for small M. As a way out, York P, |^ proposed to replace the above 
rigorous definition of the EH by an approximate criterion which is local in v and does not 
require the "teleological" information about M{v) at late v. It is supposed to be valid for 
small L. In the following we analyze 1)4.11) both using the approximate criterion and the 
exact definition of the EH. Since we have an explicit prediction for the final stages of the 
evaporation process we are in a position to determine the EH exactly, and thus to assess the 
validity and precision of York's approximate "working definition" . 

Let us begin with this working definition, applied to the RG improved spacetime. While 
according to Eq. ()4.H) . the "velocity" r{v) is strictly negative inside the TLS, the acceleration 
r{v) can have either sign there. York proposes to approximately identify r-Euiv) with the 
radius where the light ray's "acceleration" vanishes: r(EH) = 0. In pictorial terms one could 
think of this condition as separating the light rays curving downward (f < 0) from those 
curving upward (f > 0). Clearly the former (latter) are very unlikely (likely) to subsequently 




(4.1) 
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FIG. 5: Light rays of the outging null congruence with f > 0, r = 0, and r < at vi. Those with 
r > are most likely to eventually cross the TLS. 

cross the TLS, even if this is not a rigorous criterion, of course. (See Fig.(jS)) for a schematic 
sketch.) At least in the Hawking regime, photons are imprisoned by the approximate horizon 
for times long compared to the dynamical time scale of the evaporating hole 
Taking a second f-derivative of 1)4.11) we obtain the "acceleration" 

(4.2) 



L{v)^^ + G{r)M{v)^ - G'{r)M{v)- 



Here we assumed that G{r) does not have a parametric dependence on M, which is actually 
true for 7 = 0: 

G{r) = (4.3) 



r2 + uGo 

(For general 7 there appears an additional term oc dG/dM in ()4.2|) which is irrelevant 
qualitatively.) Eq. ()4.2|) tells us that, when r = 0, the radius r, the velocity r, and the time 
V are related by 

-G{r)M{v) 



-L G{r) 



- M{v)G'{r) 



(4.4) 



Since teh = '"tls if the hole would not radiate, and since higher orders in the luminosity are 
neglected, the difference ttls ~ ^eh which we would like to compute is of order L. For this 
reason we may replace r = teh on the RHS of (|4.4|) with ttls = the error being of 

order LP". 



-2LG{r^ 



2G{r+)M{v) 
r+ 



- 2M{v)G'{r^ 
15 



-2LG(r^ 



[1 - 2M{v)G'{r-+)] 



(4.5) 



In the second equality of ()4.5p we used that 2G{r+)M = r+ which follows from f{r+) = 0. 
Eq. ()4.5p provides us with the "velocity" r at the point where the acceleration vanishes. 
The corresponding coordinate, r^n, is the approximate location of the EH. We obtain it 
by using r = /(r)/2 in order to rewrite the LHS of (j4.5|) in the form /(rEH)/2, and then 
inverting the function /. This inversion is easy to perform since, again, we may expand in 
rEH -r+ = 0{L): 

fivEu) = f{r+) + (rEH - r+)/' (r+) + 0{L'). 
Hence, to order L, r(EH) = (rsH — ''^+)/'('"+)/2 which yields, together with ()4.5|) . 

_ ALG{r+) 
""^ f(r+)[l-2MG'K)] ^ 

Differentiating /(r) = 1 — 2G{r)M/r and using 2G{r^)M = r+ one obtains the relation 
f'{r-^) = [1— 2MG'(r+)]/r_|_. This leads to the following explicit formula for the f-dependence 
of Teh: 

Eq. ()4.7p is valid for an arbitrary r-dependence of G, still. If G = const we recover York's 
result Teh = r+[l- AGqL] with r+ = 2GqM. 

Let us now specialize for the function G{r) motivated by quantum gravity, Eq. (j4.3|) . 
This r-dependence of Newton's constant implies that, at the TLS, G{r^) = r^/{2M) and 
(j"(r_|_) = u3/(2M^r_|_). These equations were simplified using the relation + ujGq = 
2GoMr+ which is equivalent to = with ()4.3|) . They lead to the following result for 

the position of the EH: 

TEH = r+ [1 - 4 Go L{M) F (fi)] (4.8a) 



The radius r_|_ is given by Eq. (jl.7p or (j3.12j) . respectively, and M and fl = M^^/M{vy 
are understood to be functions of v, of course. The correction factor Y{Q) measures the 
deviations from the semi classical result; we have Y{Q) ^ 1 everywhere in the Hawking 
regime {Q ~ 0). For Q y 1 the corrections become large; in fact, Y{Q) diverges at Q = 1. 
Nevertheless, the product LY vanishes for Q 1, M \ M^, which becomes clear when we 
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FIG. 6: The EH determined numerically (solid thick line) and by York's approximate criterion 
(dashed hne) for M/Mcr = 2. The solid thin line is the AH. 



insert (fTTnil into (gSaj): 



'"EH 



2a fi(l - Vt) 



(47r)325 [1 + ^/1^] 
(47r)3 MV 



(4.9a) 



(4.9b) 



This is our final result for the radius of the event horizon, as given by York's approximate 
local criterion. In the early stages of the evaporation we recover the semiclassical result. For 
w — >■ oo however, M{v) M^r, and both r-^n and r^{v) approach the same limiting value 



asymptotically: = r±{Mcr) = VooGq = vuipi- This behavior can be seen in Fig.©. 

As we anticipated, we now see that teh is indeed smaller than r+ during the entire evap- 
oration process, the EH is inside the (outer) TLS, thus giving rise to a quantum ergosphere. 

It is an important observation that, while teh is smaller than r_|_, it is always larger 
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than r_, provided the radiation effects are small compared to the quantum gravity effects 
in an appropriate sense. Comparing the radius ()4.9|) to r_ as given by ()1.7|) we find that 
rEuiM) > r_{M) provided 

^ < (4.10) 



For < f2 < 1 the RHS of ()4.10|) is bounded below by the constant 3^/3/2. As a result, teh 
is larger than r_ during the entire evaporation process provided 

a < ^V3{47if u (4.11) 

If ()4.1H) is satisfied, we have r+{v) > r-Eu{v) > r_(t>) for any finite v, and the EH touches 
both the outer and the inner TLS only in the limit f ^ oo where r+, teh, 'T'- —>■ i^cr- 

If (14. lip is violated our method is inapplicable, most probably, and the improved Vaidya 
metric is not a reliable description of the spacetime structure. This metric is valid to first 
order in L only, which means that the dimensionless luminosity has to be much smaller than 
unity, GqL <^ 1. If so, Eq. p.lOj) implies (T/[(47r)^C;] ^ 1, and ()4.1H) is indeed satisfied. 
(Note that MlJuj"^ = Gq^/lu.) It is reassuring to see that, for pure gravity, Eq. ()4.1H) is 
fulfilled with a very broad margin; the expansion parameter (T/[(47r)^iI;] assumes the tiny 
value 8.3 ■ 10~^/(D in this case, with u = 0(1). 

Let us return to the exact definition of the EH now. We compute it for our model by 
numerically solving ()4.1|) for a set of initial conditions r(0). For a certain range of r(0)'s 
the trajectories will ultimately cross the TLS and escape to infinity, while the others remain 
at radii below ttls for arbitrarily late advanced times. The light ray separating those two 
classes of solutions defines the EH, a null hypersurface by construction. 

In Fig.(0)) we show various trajectories of either class, as well as the EH, determined 
numerically from its exact definition. We have numerically evolved several initial conditions 
r(0) in order to determine the boundary between trapped and escaping null geodesies in 
the limit v ^ oo. In Fig.® the true EH is compared to the prediction of York's criterion. 
Obviously the latter provides a rather accurate approximation to the true horizon. 

The global structure of the spacetime is depicted in the conformal diagram in Fig.Q. 
Region I is a flat spacetime, while at V = Vq {V is the Kruskal advanced time coordinate, 
defined as V = — exp{—K,v) being k the surface gravity of the outer horizon) an imploding 
null shell is present (strictly speaking it must have a negative tension in order to balance 



the flux of negative energy on its future side 
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). Region II is the evaporating black hole 



FIG. 7: The conformal diagram of the evaporating quantum black hole: region I is a flat spacetime, 
and region II is the evaporating BH spacetime, EH is the event horizon, CH is the inner (Cauchy) 
horizon, and A is the apparent horizon. 



spacetime. The AH is a timelike hypersurface which "meets" the EH at future null infinity in 
the conformal diagram. The null ray which is tangent to the earliest portion of the apparent 
horizon A would have been the EH if the hole were not radiating. The final state of the 
black hole is an extremal black hole whose inner and outer horizons have the same radius 
(r = Tcr) and are located at the event horizon EH and the inner (Cauchy) horizon CH in 



It is instructive to compare the areas A of the various horizons. They are defined 
by intersecting the EH, AH, and TLS with the incoming null surfaces v = const. Thus 
A.TLs{v) = 47rr+(t')^ and ^eh('w) = 47rrEH('w)^- From Eq. p.7j) we obtain for ^tls = ^ah 



Fig.®. 



2 



(4.12) 



and the approximate result ()4.9a|) for the event horizon implies 



-4eh — -^TLS 1 



4a VL{1 - Q) 



(4.13) 



{A-Kfu 1 + v^r^ 
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where a term of second order ma/ [AtiYuo has been neglected. The difference 5 A = .Atls — 
Aer is given by 

5A= ^ll,{l-^)[l + Vl^^ (4.14) 

During the early stages of the evaporation process, 5 A ~ o"Zpj/(27r^) = 128ttB£pi which 
coincides with the known result O for the Hawking regime, while 6 A vanishes proportional 

PR 

to (M^ — M^j.) —>■ ioT V ^ oo. It had been emphasized by York [2, y] that in the Hawking 
regime he considered, 6 A is a universal {i.e. M independent) quantity which depends only 
on (T, thus counting the degrees of freedom of the field quanta which can be evaporated off. 
Looking at Eq. ()4.14|l we see that this universality does not persist beyond the semiclassical 
approximation. 



V. CONCLUSION 



The renormalization group improvement of black hole spacetimes according to Quantum 
Einstein Gravity leads to concrete predictions on the final state of the evaporation process. 
Unlike previous studies based on ad hoc modifications of the equation of state of matter at 
very high (Planckian) densities or models based on loop quantum gravity Q, 

the mass of the remnant can be calculated explicitly: = V^ipi- Its precise value is 
determined by the value of u which is a measurable quantitv in principle j^. No naked 
singularity forms, at variance with the paradigm proposed in 4^, |^ , so that the remnant is 
a mini-black hole of Planckian size. On the other hand, it is intriguing to note that exactly 
solvable semiclassical gravity-dilaton models predict a final state described by an extremal 
configuration which is reached in an infinite amount of time [54]. (See also and js^] 



for 
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58| for an approach based on special resummations 



analogous semiclassical models, and 
of higher order graviton loops. ) 

It would be interesting to investigate the possible astrophysical implications of a pop- 
ulation of stable Planck size mini-black holes produced in the Early Universe or by the 
interaction of cosmic rays with the interstellar medium 
issues in a subsequent publication. 
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0- 



We hope to address these 



20 



VI. ACKNOWLEDGMENTS 



We would like to thank Werner Israel for helpful suggestions on an earher version of this 
manuscript. A.B. would also hke to thank the Department of Physics of Mainz University 
and M.R. the Catania Astrophysical Observatory for the financial support and for the cordial 
hospitality extended to them while this work was in progress. We are also grateful to INFN, 
Sezione di Catania for financial support. 



[1] S. Hawking, Commun. Math. Phys. 43, 199 (1975). 

[2] J. York, in Quantvm Theory of Gravity: Essays in Honor of the Sixtieth Birthday of 

Bryce S. DeWitt, edited by S. Christensen (Adam Hilger, Bristol, 1984). 
[3] J. York, Phys. Rev. D 28, 2929 (1983). 
[4] J. York, Phys. Rev. D 31, 775 (1985). 
[5] P. Hajicek and W. Israel, Phys. Lett. A80, 9 (1980). 
[6] E. Poisson, A Relativist's Toolkit (Cambridge University Press, 2004). 
[7] M. Renter, Phys. Rev. D 57, 971 (1998). 

[8] D. Don and R. Percacci, Class. Quantum Grav. 15, 3449 (1998). 

[9] O. Lauscher and M. Renter, Phys. Rev. D 65, 025013 (2002). 
[10] M. Renter and F. Saneressig, Phys. Rev. D 65, 065016 (2002). 
[11] O. Lauscher and M. Renter, Class. Quantum Grav. 19, 483 (2002). 
[12] W. Souma, Prog. Theor. Phys. 102, 181 (1999). 
[13] R. Percacci and D. Perini, Phys. Rev. D 68, 044018 (2003). 
[14] M. Reuter and F. Saueressig, Phys. Rev. D 66, 125001 (2002). 
[15] D. F. Litim, Phys. Rev. Lett. 92, 201301 (2004). 
[16] M. Niedermaicr, Nucl. Phys. B673, 131 (2003). 
[17] O. Lauscher and M. Reuter, Phys. Rev. D 66, 025026 (2002). 
[18] R. Percacci and D. Perini, Phys. Rev. D 67, 081503 (2003). 
[19] M. Reuter and F. Saueressig, Fortsch. Phys. 52, 650 (2004). 
[20] M. Niedermaier, JHEP 12, 066 (2002). 
[21] P. Forgacs and M. Niedermaier (2002), hep-th/0207028. 



21 



[22] O. Lauscher and M. Reuter, Int. J. Mod. Phys. A 17, 993 (2002). 

[23] O. Lauscher and M. Reuter, JHEP 10, 050 (2005). 

[24] O. Lauscher and M. Reuter (2005), hep-th/0511260. 

[25] M. Reuter and J.-M. Schwindt, JHEP 0601, 070 (2006). 

[26] A. Bonanno and M. Reuter, JHEP 0502, 035 (2005). 

[27] A. Bonanno and M. Reuter, Phys. Rev. D 60, 084011 (1999). 

[28] A. Bonanno and M. Reuter, Phys. Rev. D 62, 043008 (2000). 

[29] A. Bonanno and M. Renter, Phys. Rev. D 65, 043508 (2002). 

[30] A. Bonanno and M. Reuter, Phys. Lett. B527, 9 (2002). 

[31] E. Bentivegna, A. Bonanno, and M. Reuter, JCAP 0401, 001 (2004). 

[32] A. Bonanno, G. Esposito, and C. Rubano, Gen. Rel. Grav. 35, 1899 (2003). 

[33] T. Tsuneyama (2004), hep-th/0401110. 

[34] M. Reuter and H. Weyer, Phys. Rev. D 69, 104022 (2004). 

[35] M. Reuter and H. Weyer, JCAP 12, 1 (2004). 

[36] M. Reuter and H. Weyer, Phys. Rev. D 70, 124028 (2004). 

[37] F. Bauer, Class. Quantum Grav. 22, 3533 (2005). 

[38] C. Wetterich, Phys. Lett. B301, 90 (1993). 

[39] M. Reuter and C. Wetterich, Nucl. Phys. B417, 181 (1994). 

[40] M. Reuter and C. Wetterich, Nucl. Phys. B427, 291 (1994). 

[41] M. Reuter and C. Wetterich, Nucl. Phys. B391, 147 (1993). 

[42] M. Reuter and C. Wetterich, Nucl. Phys. B408, 91 (1993). 

[43] J. Bardeen, Phys. Rev. Lett. 23, 2823 (1981). 

[44] W. Hiscock, Phys. Rev. D 23, 2813 (1981). 

[45] W. Hiscock, Phys. Rev. D 23, 2823 (1981). 

[46] N. Birrell and P. Davies, Quantum Fields in Curved Space (Cambridge University Press, 1982). 

[47] S. Hawking and G. Ellis, The Large Scale Structure of Space-Time (Cambridge University 
Press, 1973). 

[48] C. Barrabes and W. Israel, Phys. Rev. D 43, 1129 (1991). 

[49] V. P. Prolov, M. A. Markov, and V. F. Mukhanov, Phys. Rev. D41, 383 (1990). 

[50] I. Dymnikova, Int. J. Mod. Phys. D12, 1015 (2003). 

[51] S. A. Hayward (2005), gr-qc/0506126. 

22 



[52] A. Ashtekar and M. Bojowald, Class. Quant. Grav. 22, 3349 (2005). 
[53] J. F. Donoghue, Phys. Rev. Lett. 72, 2996 (1994). 

[54] A. Fabbri, D. J. Navarro, and J. Navarro-Salas, Phys. Rev. Lett. 85, 2434 (2000). 

[55] D. Grumiller, JCAP 0405, 005 (2004). 

[56] T. Jacobson, Phys. Rev. D 57, 4890 (1998). 

[57] B. F. L. Ward, Int. J. Mod. Phys. A20, 3502 (2005). 

[58] B. F. L. Ward (2005), hep-ph/0502104. 

[59] A. Barrau, C. Feron, and J. Grain, Astrophys. J. 630, 1015 (2005). 

[60] A. Barrau, D. Blais, G. Boudoul, and D. Polarski, Phys. Lett. B551, 218 (2003). 



23 



